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Abstract 

We investigate the p(^, K + )A and p{^,K + )TP reactions within a coupled-channels effective- 
ly'! 

Lagrangian method which is based on the K-matrix approach. The two-body final channels 
Q-i' included are nN, nN, (j)N, pN, jN, KA, and KT,. Non-resonant meson-baryon interactions 



43 



are included in the model via nucleon intermediate states in the s- and u-channels and meson 



£NJ ■ exchanges in the i-channel amplitude and the u-channel resonances. The nucleon resonances 

> 

5n(1535), 5n(1650), 5 3 i(1620), Pu(1440), Pn(1710), P 13 (1720), P 33 (1232), P 33 (1600), D 13 (1520), 
Z?i 3 (1700), and D 33 (1700) are included explicitly in the calculations. With a single parameter set 
which was derived earlier from our analysis of the n meson photoproduction, the model describes 
ON . well all the available cross section and polarization data of the SAPHIR collaboration for the two 

investigated channels. The description of the data of the CLAS collaboration, however, is not of 



^ ■ the same quality. In contrast to some previous studies, we do not find any compelling need for in- 

eluding a I?i 3 state with mass of around 2.0 GeV in order to reproduce the data for the ^(7, K + )A 
reaction at photon energies corresponding to the invariant mass around 1.9 GeV. 

PACS numbers: 13.60.ie, 13. 75. Cs, 11.80. - m, 12.40. Vv 
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I. INTRODUCTION 



Recently, several laboratories have devoted considerable amount of effort to the inves- 
tigation of the photoproduction of strangeness on the nucleon. Experiments performed at 
JLab-CLAS HQ, ELSA-SAPHIR 3, 4] and SPring8/LEPS ^ 6] have produced high quality 
data on the associate strangeness production reactions p(7, K + )A and ^(7, K + )TP covering 
the photon energy regime from threshold to upto 3.0 GeV. Furthermore, data on single and 
double polarization |7| observables have also become available. The motivation behind these 
studies has been the fact that a considerable part of the excitation spectrum of the nucleon 
can, in principle, participate in the production process of associated strangeness (K + A and 
K + S°) through these reactions because even at the threshold they involve invariant masses 
that exceed those of several baryonic resonances. It is hoped that these studies could help 
in investigating the so-called missing baryonic resonances that are predicted by the quark 
models (see, e.g., Ref. js]) but not observed in non-strange meson photoproduction. Some 
of these resonances may couple strongly to the KA and KT, channels. 

The determination of the properties of the nucleon resonances (e.g., their masses, widths, 
and coupling constants to various decay channels) is an important issue in hadron physics. 
This will help in testing the predictions of lattice quantum chromodynamics (LQCD) which 
is the only theory that tries to calculate these properties from first principles. Even though, 
the requirement of computational power is enormous for their numerical realization, such 
calculations have started to provide results for the properties of nucleon ground as well 



as excited states 



.0 



13j. Furthermore, reliable nucleon resonance data are important for 



testing the "quantum chromodynamics (QCD) based" quark models of the nucleon (see, 
e.g., 14)) and also the dynamical coupled- channels models of baryonic resonances 15]. 

One of the major challenges of this field is the extraction of reliable information about 
nucleon resonance properties from the photoproduction data. In experiments, these reso- 
nances are excited as intermediate states before decaying into the final meson and baryon 
channels. A good description of intermediate-energy scattering is still not amenable to the 
LQCD calculations. Therefore, at this stage the effective methods are usually employed 
to describe the dynamics of the meson production reactions. Baryon resonance states are 
included explicitely in these approaches and their properties are obtained by comparing the 
predictions of the theory with the experimental data 161—121 1 . 
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In order to determine resonance properties reliably from the experimental measurements 
one requires a model that can analyze the different reactions over the entire energy range 
using a single Lagrangian density that generates all non-resonance contributions from Born, 
u- and t-channel contributions without introducing new parameters. At the same time, the 
Lagrangian should also satisfy the symmetries of the fundamental theory (i.e. QCD) while 
retaining only mesons and baryons as effective degrees of freedom. 

The coupled-channels method within the K-matrix approach 22H28| provides a way to 
analyze simultaneously all reaction data for a multitude of observables in different reaction 
channels while respecting the constraints described above. This method is attractive because 
it is based on an effective-Lagrangian framework that is gauge invariant and is consistent 
with chiral symmetry. It also provides a convenient way of imposing the unitarity constraint. 
This results from the Bethe-Saltpeter equation in the approximation where particles forming 
the loop are taken on the mass shell i.e. only the discontinuity part of the loop integral is 
retained. The S matrix in this approach is unitary provided the K-matrix is taken to be 
real and Hermitian. 

Alternatively, the dynamical coupled channels models within the Hamiltonian formalism 



have also been used to describe the meson-production reactions 



29-35 



such as Kaon-Maid 



361 ] and Saclay-Lyon [37J have been utilized in Ref. 



Isobaric models 



38 



391 to describe 



the 0(7, K + )A reaction. This reaction has also been studied within a variety of tree-level 

nn nn n 

isobar models [4CH43I] and in the quark models j44j-|46j . In Ref. |47| a gauge- invariant chiral 
unitary framework is used to describe the photoproduction data of both SAPHIR and CLAS 
collaborations. 

The main objective of this paper is to study photoproduction reactions ^(7, K + )A and 
p(7, K + )TP for photon energies ranging from threshold to about 3 GeV in a coupled-channels 



formalism of Refs. 



25 



27 



which is based on the K-matrix approach. This is an effective 
Lagrangian model which is gauge invariant and obeys the low-energy theorem. We aim 
at describing simultaneously the data on total and differential cross sections as well as on 
polarization observables for both the reactions within the same framework with a single 
parameter set. The A and S° hyperons have isospins of and 1, respectively. Thus, the 
intermediate states of K + A have the isospin | only (A*) while those of K + TP can have both 
\ and I isospins (A* and A). Therefore, a combined description of all the available data for 
both these channels is indeed quite interesting and is a challenge to any theoretical model. 
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We would like to add that a subset of SAPHIR and CLAS data for the p(7, K + )A reaction 



has been investigated previously in Ref. 



261 ] within the Giessen model which is also a coupled- 



channels effective-Lagrangian K-matrix approach. Despite some differences in details (see 



the discussions in Ref. 



28]), our method is similar to that of the Giessen model. The 



calculations presented in Ref. 26] have reproduced the p(-j, K + )A data well with a slight 
preference to the SAPHIR data. In the present work we have attempted to describe the 
p('~f,K + )Tj reaction also along with the p( r y,K + )h one within the same framework using 
the same parameter set which already makes our analysis self-contained. Nevertheless, by 
comparing our results with those of the Giessen model we expect to gain further insight in 
the mechanism of the strangeness photoproduction off proton. 

Our paper is organized in the following way. An overview of our model is given in section 
II. This consists of a short discussion of the K-matrix formalism, the model space and the 
channels included, the Lagrangians, and the form factors. Our results and their discussions 
are presented in section III. Summary and conclusions of our work are presented in section 
IV. Finally form of Lagrangians at various vertices are given in appendix A. 



II. DESCRIPTION OF THE MODEL 



This work is based on an effective-Lagrangian model. The kernel in the K-matrix approach 
is built by using the effective Lagrangian as given in Appendix A. We have taken into 
account contributions from (i) the nucleon Born term, (ii) t-channel exchanges of mesons, 
(iii) nucleon and resonance terms in the w-channel, and (iv) baryonic resonance in the s- 
channel (see Fig. 1). The sum of amplitudes (i), (ii) and (iii) is termed as the background 
contribution. As is discussed below, this approach allows to account for coupled-channels 
effects while preserving many symmetries of a full field-theoretical method. 



A. K-matrix model 



The coupled-channels (or re-scattering) effects are included in our model via the K-matrix 
formalism. In this section we present a short overview of this approach; a more detailed 



description can be found in Refs. 



23 



25 



48 



49]. 
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FIG. 1: [color online] Feynman diagrams included in this work. First row: s- and w-channel 
diagrams with propagating final state baryons (B = N, A, S) or intermediate state resonances 
(A, N*). M stands for mesons included in the model space. Second row: i-channel contributions 
with propagating asymptotic and intermediate mesons, and the contact term required by the gauge 
invar iance. 



In the K-matrix formalism the scattering matrix is written as 

It is easy to check that the resulting scattering amplitude 5=1 + 2iT is unitary provided 
that K is Hermitian. The construction in Eq. (CQ) can be regarded as the re-summation of 
an infinite series of loop diagrams by making a series expansion, 

T = K + iKK + z 2 KKK + ■ ■ • . (2) 

The product of two K-matrices can be rewritten as a sum of different one-loop contributions 
(three- and four-point vertex and self-energy corrections) depending on Feynman diagrams 
that are included in the kernel K. However, not the entire spectrum of loop corrections 
present in a true field-theoretical approach, is generated in this way and the missing ones 
should be accounted for in the kernel. In constructing the kernel, care should be taken to 
avoid double counting. For this reason we include in the kernel tree-level diagrams only 
[Figs. l(a)-l(c)], modified with form-factors and contact terms [Fig. 1(d)]. The contact 



TABLE I: Baryon states included in the calculation of the kernel with their coupling constants. 
The column labeled WD lists the decay width to states outside the model space. The columns 
labeled M and WD are in units of GeV. See text for a discussion on the signs of the coupling 



constants. 



Lij 


M 


WD 


gNn 


<?P7 


2 


9K A 


9KY. 


3JVtj 


5n(1535) 


1.525 


0.0 


0.6 


-0.60 




0.1 


0.0 


2.2 


Sn(1650) 


1.690 


0.030 


1.0 


-0.45 




-0.1 


0.0 


-0.8 


531(1620) 


1.630 


0.100 


3.7 


-0.12 






-0.8 




Pii(1440) 


1.520 


0.200 


5.5 


0.65 




0.0 


-2.0 


0.0 


Fu(1710) 


1.850 


0.300 


3.0 


0.25 




0.0 


-3.0 


2.0 


Fl 3 (1720) 


1.750 


0.300 


0.12 


-0.75 


0.25 


-0.05 


0.0 


0.12 


^33(1230) 


1.230 


0.0 


1.7 


-2.2 


-2.7 




0.0 




^33(1600) 


1.855 


0.150 


0.0 


-0.4 


-0.6 




0.55 




£>i 3 (1520) 


1.515 


0.050 


1.2 


2.6 


2.5 


2.0 


0.0 


1.2 


£>i 3 (1700) 


1.700 


0.090 


0.0 


-0.5 


0.0 


0.0 


0.3 


-0.04 


^33(1700) 


1.670 


0.250 


0.8 


1.5 


0.6 




-3.0 





terms (or four-point vertices) ensure gauge invariance of the model and express model- 
dependence in working with form factors (see Section III CI) . Inclusion of both s- and u- 
channel diagrams [Figs 1(a) and 1(b), respectively] in the kernel insures the compliance 
with crossing symmetry. 

To be more specific, the loop corrections generated in the K-matrix procedure include only 



50 



This is 



those diagrams which correspond to two on-mass-shell particles in the loop 
the minimal set of diagrams one has to include to ensure two-particle unitarity. Thus, all 
diagrams that are not two particle reducible (e.g. 7iV — > 2nN channels), are not included. 
In addition, only the convergent pole contributions i.e. the imaginary parts of the loop 
correction, are generated. The omitted real parts are important to guarantee analyticity 
of the amplitude and may have complicated cusp-like structures at energies where other 
reaction channels open. In principle, these can be included as form factors as is done in the 
dressed K-matrix procedure 50j, [52(]. We have chosen to work with purely phenomenological 



form factors for the reason of simplicity. An alternative procedure to account for the real- 
loop corrections is offered in Refs. 53m5j] . 

The strength of the K-matrix procedure is that in spite of its simplicity, several symmetries 



are obeyed by it 



481 ] . As was already noted, the resulting amplitude is unitary provided 



that K is Hermitian, and it obeys gauge invariance provided the kernel is gauge- invariant. 
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In addition, the scattering amplitude complies with crossing symmetry when the kernel is 
crossing symmetric. This property is crucial for a proper behavior of the scattering amplitude 



in the low-energy limit 5l|, |56[. Coupled-channels effects are automatically accounted for 



by this approach for the channels explicitly included into the K-matrix as the final states. 

As a result of this channel coupling, the resonances generate widths which are compatible 
with their decays to channels included in the model space. For some resonances, such as the 
A and the 5*11(1535), this corresponds to their total width. Other resonances, particularly 
the high lying ones, may have important decay branches to states that are not included in 
the model basis. To account for this in our calculations, we have added an explicit dissipative 
part to the corresponding propagators. The magnitudes of these widths are equivalent to 
decay widths of the resonances to states outside of our model space. 

The resonances which are taken into account in building the kernel are summarized in 
Table HI In the current work we limit ourselves to the spin-| and spin-| resonances as in 
this energy regime higher spin resonances are known 26[ to give only a minor contribution 
to the K + A channel. Spin-| resonances are included with so-called gauge- invariant vertices 
which have the property that the coupling to spin-| pieces in the Rarita-Schwinger prop- 



agator vanish 



50 



57 



5N 



We have chosen this prescription since it reduces the number 



of parameters as we do not have to deal with the off-shell couplings. The effects of these 



5s 



which we prefer, certainly within 



off-shell couplings can be absorbed in contact terms 
the context of the present work. 

The masses of the resonances given in Table [I] are bare masses and they thus may deviate 
from the values given by the Particle Data Group 59[. Higher-order effects in the K-matrix 
formalism do give rise to a (small) shift of the pole-position with respect to the bare masses. 
The masses of very broad resonances, in particular the Pn, are not well determined - values 
lying in a broad range (typically a spread of the order of a quarter of the width) give 
comparable results. The width quoted in Table H] corresponds to the partial width for decay 
to states outside our model space. The parameters as quoted in Table [Dare mostly unchanged 



as compared to those presented in previous calculations within this model 25 



27 



. The 



t-channel contributions which are included in the kernel, are summarized in Table [TTJ 
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TABLE II: Mass, spin, parity and isospin of the mesons which are included in the model. The 
rightmost column specifies in which reaction channels their f-channel contribution are taken into 
account. 



Meson 


M [GeV] 




I 


t-ch contributions 


7T 


0.135 


o - 


1 


f'v/V - 


-V riiftH ( 7r/V 


— ^ />iV ) 


K 


0.494 


o - 


1 

2 


V 7 iv 








1.019 


1- 











7] 


0.547 


o - 











P 


0.770 


1- 


1 


( 7 A -> irN), 


(7 A -> 77A), 


(AA AS), 










(AS - 


■> AS), (Ayr - 


^ AA), 










(Ntt - 


-+ At/), (Avr - 


-> Avr) 


UJ 


0.781 


1- 





(iV 7 - 


-> Avr), (A7 - 


-> At?) 


a 


0.760 


0+ 





(iV 7 - 


-> A<£), (Avr - 


-> Avr) 


K* 


0.892 


1- 


1 

2 


(A 7 - 


+ AA), (A 7 - 


4 AS), 










(AA - 


-4 At?), (AS - 


-4 Ar/), 












(Avr -> AS) 





B. Model space, channels included 

To keep the model manageable and relatively simple, we consider only stable particles 
or narrow resonances in two-body final states. The AK, EAT, N<f>, Nrj and N'-f are the 
final states of primary interest, and the Nit final state is included for its strong coupling to 
most of the resonances. Three-body final states, such as 2ttN, are not included explicitly 
as stated above. Their influence on widths of resonances is taken into account by assigning 



an additional (energy dependent) width to them 23|. To investigate the effects of coupling 



to more complicated states, we have also included the Np final state. As was shown in 



Ref. [25|, inclusion of the p channel has a strong influence on the pion sector but only a 
relatively minor effect on A and S photoproduction. 

The components of the kernel which couple the different non-electromagnetic channels 
are taken as the sum of tree-level diagrams, similar to what is used for the photon channels. 
For these other channels no additional parameters were introduced and they thus need no 
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further discussion. 



C. Form- factors &: gauge restoration 

Without the introduction of form factors, calculations with Born terms strongly overes- 
timates the cross section at higher energies. Although inclusion of coupled-channels effects 
reduces the cross section at high energies yet disagreement with the experimental data still 
persists. Therefore, the Born contribution will have to be quenched with form factors. There 
are two physical motivations for introducing form factors (or vertex functions). First of all, 
at high photon energies one may expect to become sensitive to the short-range quark struc- 
ture of the nucleon. Because this physics is not included explicitly in our model, we can 
only account for it through the introduction of phenomenological vertex functions. The sec- 
ond reason has to do with the intermediate-range effects because of meson-loop corrections 
which are not generated through the K-matrix formalism. Examples of these are given in 
Refs. 



50 



52]. 



In our approach as well as that of Ref. 



24|, the form- factors are not known a priori 



and thus they introduce certain arbitrariness in the model. In the current paper we limit 
ourselves to dipole form-factors in s-, u-, and t-channels because of their simplicity, 

A 2 

X z + {s — m z ) z 

where m is the mass of the propagating particle and A is the cut-off parameter. For ease of 
notation we introduce the subtracted form factors 

L(s)= 1 ~ Fm [ S \ (4) 
s — m z 

where F m (s) is normalized to unity on the mass-shell, F m (m 2 ) = 1, and / m (m 2 ) is finite. 

However, only in the kaon sector we use a different functional form for the it-channel 
form-factors 

Hm ^ = f\2,f U 2YA — 2" ( 5 ) 

(A z + (it — m z ) z jm z 

The argumentation for this different choice is presented in the discussion of the E 
photoproduction results in Ref. 



251 ] . Often a different functional form and cut-off values 
are introduced for the t-channel form factors. Although this can easily be motivated, it 
introduces additional model dependence and increases the number of free parameters. To 

9 



limit the overall number of parameters we have taken the same cut-off value [A = 1.2 GeV 2 , 
see Eq. (j3J)] for all form-factors except for the Born contributions in kaon channels where we 
used A = 1.0 GeV 2 . 

Inclusion of form-factors will in general break electromagnetic gauge-invariance of the 



25|, the impli- 



model. Therefore, a gauge-restoration procedure should be applied. In Ref. 
cations of various gauge-restoration procedures was studied for the 7^ — > KT, amplitude. It 
was observed that the gauge-invariance restoration procedure is model dependent which may 
give rise to strongly different Born contributions to the amplitude. Therefore, the choice 
of a procedure to be adopted is guided by its ability to describe the experimental data. It 
was found that the gauge-restoration procedure of Davidson and Workman 60( provided 
the best description of the data on the KH photoproduction. We have used this procedure 
in the present work also. 

We note that fitting the pion-scattering and pion-photoproduction amplitudes fixes 
masses as well as pion- and photon-coupling constants for most of the resonances. This 
limits strongly the number of free parameters for the kaon-production channels. 



III. RESULTS AND DISCUSSIONS 

Our aim in this paper is to use the data base on p(~f, K + )A and p(7, K + )T,° reactions of 
both SAPHIR and CLAS collaborations to check various ingredients of our unitary coupled- 
channels field theoretic model of meson production in photon induced reactions on nucleons. 
In particular we are interested in checking to what extent a simultaneous fit to data for a 
multitude of observables for both the reactions can be obtained with a single set of input 
parameters. This is expected to provide a strong constraint on the model parameters thus 
reducing the model dependence to a minimum. It is also likely to highlight the role of 
channel couplings in various regions of photon energies because several calculations of the 
associated kaon production reactions have neglected these effects. 

We emphasize however, that even the large experimental data base may not allow to fix 
the extracted parameters uniquely within the unitary coupled-channels effective-Lagrangian 



model 24|. This is due to the fact that it is necessary to include empirical form factors 
in the model to regularize the amplitudes at higher energies. As mentioned before These 
form factors require a gauge-invariance restoration procedure which involves ambiguities. 
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FIG. 2: [color online] Comparison of the calculated total cross sections for the !_p(7, K + )K and 
p(7, K + )TP reactions with the corresponding experimental data taken from Ref. [3(. 

Nevertheless, confronting the model with a large data base in several reaction channels is 
expected to provide a means to overcome this problem. 



The parameters in the model have been adjusted 
partial wave amplitudes of Arndt et al. 61] . In Ref. 



27] to reproduce the Virginia Tech 
we have presented a comparison 



of our calculated S-, P-, and D-wave amplitudes for pion-nucleon scattering for isospins / 
= 1/2 and 3/2 channels with those of the FA08 single-energy partial wave amplitudes of 
Ref. The corresponding results for pion photoproduction and the Compton scattering 
are given in Ref. . There we noted that both real and imaginary parts of the pion-nucleon 
scattering amplitudes are described well although some differences start to show up at the 
upper limit of the energy range considered. 



The data for p( / y,K + )A and p(7, K + )H° reactions consist of total and differentia 



sections and hyperon polarizations measured at CB-ELSA (SAPHIR collaboration J3]) and 
at JLab (CLAS collaboration jl, 2|) for photon energies ranging from respective thresholds 
to about 3 GeV. Moreover, beam asymmetry data are available for 9 photon energy beans 
between 1.5 GeV to 2.4 GeV from SPring8/LEPS group. These data, therefore, cover 
not only the entire resonance region but also the region where the background contributions 
are expected to be dominant. 

In Fig. 2, we compare the results of our calculations for the total cross sections with the 
corresponding data of SAPHIR collaboration js] for p('~f,K + )A and p(7,i^ + )S° reactions. 
Photon energies (J5 7 ) range from threshold to about 2.6 GeV. The experimental cross sections 



cross 



11 



for the p(7, K + )A reaction show a steep rise as P 7 increases from threshold to about 1.1 GeV. 
The latter corresponds to a 7£> channel total invariant mass (W) of ~ 1.7 GeV which coincides 
with the masses of Sii(1650), Pn(1710) and Pi3(1720) resonances. The decrease in the cross 
sections just before the threshold for the K + Tf channel indicates a cusp due to the opening 
of this channel which has already been indicated in Refs. 25|, |62|. The data in this region 



are well reproduced by our calculations. There is also a second peak in the data at E 1 ~ 
1.5 GeV ( W ~ 1.9 GeV). Our calculations describe the data well also in this region. 



In contrast to studies reported in Ref. 



21 



38 



39 



631 ]. we do not require an additional 



baryonic resonance Pi3(1895) in order to explain the data in the second peak region. Looking 
at the cross sections for the p( , y,K + )'E reaction (right panel) one notices that the second 
maximum in the p(7, K + )A data is centered around the same value of W where the only 
peak is observed in the K + TP channel. Hence, peaks in the p(7, K + )A total cross section 
are more likely to be the consequences of unitarity and multi-channel dynamics. This is 
further supported by the coupled channel analysis of this reaction as reported in Ref. {26 1 
where the second peak in the K + A total cross section data is explained as resulting from the 
interference between background and resonance contributions and not due to the presence of 
a i?i3(1895) resonant state. The same conclusion was arrived also in the previous coupled- 



channels studies [24 1 . 

On the other hand, the total cross sections of the ^(7, K + )TP reaction rise smoothly from 
threshold to its peak at P 7 « 1.45 GeV (W ~ 1.9 GeV). The cross sections drop smoothly 
for Ej larger than this value. Our calculations are able to reproduce the data well in the 
entire region of photon energies with the exception of some far lower and some higher photon 
energies where the data are somewhat underestimated. 

In Fig. 3, we show the contribution of various resonances and background terms to the 
total cross sections of ^(7, K + )A and p(7,i^ + )E° reactions as a function of E 1 . It is clear 
from this figure that while the background contributions dominate the ^(7, K + )A cross 
sections in the entire range of photon energies, they do so only for P 7 > 1.5 GeV in case of 
the p(7, K + )TP reaction. It is interesting to note that contributions of Pi3(1720), Sii(1535), 
S n (1650), P 13 (1720) and P> 13 (1520) (depicted as P13, 511-1, 511-2, P13-1, respectively in 
the left panel of Fig. 3) resonances peak at about the same value of P 7 (~ 1.1 GeV) in the 
p(7, K + )A total cross section. In the region around this energy, the resonance contributions 
are comparable to those of the background terms and they combine together to constitute 
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FIG. 3: [color online] Partial wave decomposition of the calculated total cross sections for the same 
reactions as those shown in Fig. 2. Contributions of different resonances are shown by various 
curves as indicated in the figure. Also shown are the background contributions which consist of 
Born and u- and t-channel terms. Experimental data are taken from Ref. 3|. 



the structure of the first peak in the data. Furthermore, no one resonance individually 



dominates in this region which is in contrast to the results of Ref. 



26[. The contributions 



of the Pn(1710) resonance are very weak and are not included in our study. This is in 



agreement with the results Ref. 



26] . We have also not included a third Su resonance with 



mass and width around 1.780 GeV and 0.28 GeV, respecti vely which was considered in 
descriptions of p(7, rj)p and p(7, K + )A reactions in Refs. (33), 46| . 

The peak region in the p(7,i^ + )E° total cross section, on the other hand, is dominated 
by the contributions from the spin-|, isospin-| P 33 (1600) resonance. Apart from the back- 
ground and to a lesser extent .033(1 700) terms other resonances are almost unimportant 
in this region. Furthermore, magnitudes of Di 3 (1700), Pn(1440) and Pn(1710) resonances 
(depicted as P13-2, PI 1-1 and PI 1-2, respectively in the right panel of Fig. 3) are compar- 
atively small in the entire range of photon energies. However, for E 1 very close to threshold 
the S31 (1620) resonance is most important. 
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FIG. 4: [color online] Comparison of the calculated and experimental differential cross sections 
for the p(7, K + )A and p(^, reactions as a function of the cosine of the K + cm. angle for 

photon energies 0.94 GeV < E 1 < 1.4 GeV and and 1.05 GeV < E 1 < 1.55 GeV, respectively. The 
energy bin is indicated in each graph in GeV. The experimental data are taken from Ref. 3(]. 



For both the reactions shown in Fig. 3, we note that the total cross sections beyond 2 GeV 
are almost solely governed by the contributions of the background terms. In this region all 
resonance contributions are small and comparable to each other. 

Differential cross sections (DCS) provide more valuable information about the reaction 
mechanism. They reflect the quantum number of the excited state (baryonic resonance) 
when the cross section is dominated by it. DCS include terms that weigh the interference 
terms of various components of the amplitude with the outgoing K + angles. Therefore, 
the structure of interference terms could highlight the contributions of different resonances 
in different angular regions. For ^(7, K + )A and p(7,i^ + )S° reactions DCS data of ELSA- 
SAPHIR group exist for 36 and 35 photon energy bins, respectively in the range of respective 
thresholds to about 2.60 GeV covering a wide range of K + center of mass (cm.) angles jsj. 
In the left panels of Figs. 4-6 we show comparisons of our calculations for DCS with the 
corresponding SAPHIR data for the ^(7, K + )A reaction for energy bins in the range of 
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FIG. 5: [color online] Same as that shown in Fig. 4 but for photon energies 1.4 GeV < E~ < 2.0 
GeV and 1.55 GeV < £7 7 < 2.15 GeV, respectively. The energy bin is indicated in each graph in 
GeV. The experimental data are taken from Ref. |3j|. 
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0.9 GeV - 1.4 GeV, 1.4 GeV - 2.0 GeV, and 2.0 GeV - 2.6 GeV, respectively, while in the 
right panels the same are shown for the p(^y, K + )TP reaction for photon energy bins in the 
range of 1.05 GeV - 1.55 GeV, 1.55 GeV - 2.15 GeV, and 2.15 GeV - 2.55 GeV, respectively. 



It is seen that the differential cross sections are flat as a function of kaon angle near the 
respective thresholds which signifies the dominance of 5*-wave resonances near these energies. 
As E~f rises further the DCS develop a significant forward peaking which is consistent with 
the domination of the background terms or the interference between background and s- 
channel resonance contributions. At still higher energies, the data show a tendency of a 
slow rise at the extreme forward angles. 

It is clear that our model describes general trends of the data well in the complete photon 
energy regime of the SAPHIR measurement. However, a few specific details of the data are 
missed for some energy bins. In the regions of 1.0-1.20 GeV and 1.45-1.8 GeV, the angular 
distributions at the extreme forward angles are not properly reproduced for the p(j, K + )A 
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FIG. 6: [color online] Same as that shown in Fig. 4 but for photon energies 2.0 GeV < E~ < 2.6 
GeV and 2.15 GeV < E^ < 2.55 GeV, respectively. The energy bin is indicated in each graph in 
GeV. The experimental data are taken from Ref. [3(. 

reaction. However, those of the p( , ~f,K + )'E reaction are well reproduced in these regions. 
At photon energies > 2.0 GeV, the p(7, K + )A cross sections show a tendency of peaking at 
extreme forward angles while the p(7, K + )TP data do not seem to do so. Within statistical 
errors our calculations are consistent with this trend of the data, although for the K + A 
channel the agreement is of a lesser quality as compared to that for the K + TP one. In view 
of the fact that our background terms include both K and K* exchange diagrams with the 
same couplings for both the reactions, we have obtained a reasonably good agreement with 
the data. 

There are some discrepancies between the data of ELSA-SAPHIR and CLAS collabora- 
tions. The CLAS group has reported consistently larger cross sections at most kaon angles 
for E 1 > 1.19 GeV. For certain forward angles and photon energy bins one notices a large 
difference in the SAPHIR and CLAS data for DCS. Therefore, a simultaneous description of 
the data of two collaborations has often been problematic for the theoretical models 
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Our model is no exception to this. We show in Figs. 7 and 8, comparisons between our cal- 
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FIG. 7: [color online] Comparison of the calculated and experimental differential cross sections 
from the CLAS collaboration jl] (shown by solid blue squares) for the p(j, K + )K and p{^,K + )TP 
reactions as a function of the cosine of the K + cm. angle for selected bins of photon energy < 1.5. 
The data of the SAPHIR collaboration (taken from Ref. [J]) at nearby energies are also shown for 
comparison (by solid black circles). 

culations and the CLAS data for differential cross sections for p(7,i^ + )A and p(7,i^ + )E° 
reactions for 12 chosen photon energy bins. In these figures the SAPHIR data points are also 
given for the comparison purpose. Since, the CLAS data are given in energy and angular 
bins different from those of the SAPHIR one, we have chosen those bins which are nearly 
equal in two cases. 

We note that CLAS and SAPHIR data are identical for E 1 < 1.1 GeV for both reac- 
tions. However, at larger energies considerable difference are seen between the two data 
sets particularly at forward angles. The CLAS DCS data for the K + A channel are more 
forward peaked for E 1 & 1.2 - 1.4 GeV as compared to the SAPHIR data as well as our 
calculated cross sections. For the K + TP case, the CLAS DCS are consistently larger than 
the SAPHIR ones for all angles except for the backward ones for between 1.25 to 1.5 
GeV. Our calculations are unable to reproduce this feature. 

For E 1 between 1.4 -1.7 the K + A channel CLAS data have the tendency of backward 
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FIG. 8: [color online] Same as that shown in Fig. 7 but for selected bins of photon energy > 1.5 
GeV. 
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peaking as well along with the stronger forward peaking which is again in contrast to the 
SAPHIR data as well as our calculations. The backward peaking of the CLAS data dis- 
appears for photon energies between 2.0 - 2.6 GeV, however, stronger forward peaking still 
remains and our calculations are unable to reproduce it fully. For the p( r f,K + )YP reaction 
the CLAS data show some backward peaking at photon energies between 2.2 - 2.4 GeV. At 
other energies the differences between the data of the two groups are less noticeable for this 
channel. Our calculations are in better agreement with the CLAS data for these cases. 

Resonance structure in the s-channel should appear more clearly in the W dependence 
of the differential cross sections at various K + angles. In Fig. 9 (Fig. 10) we show the DCS 
for p(7, K + )A and p(7, K + )TP reactions as a function of W for positive (negative) values of 
COS(#k+). The experimental data are from the SAPHIR collaboration. From these figures 
we note that overall shapes of the W distributions are reproduced well by our calculations 
for all the angles for both the reactions. 

Nevertheless, we also notice that for the case of the K + A production our calculations un- 
derestimate the data somewhat for W between 1.7-2.0 GeV at forward angles [for COS(9k+) 
= 0.75, 0.85, and 0.95]. As bulk of the forward peaking is due to t-channel exchanges which 



18 



0.8 
0.6 
0.4 
0.2 


0.6 
0.4 
0.2 


0.6 
0.4 
0.2 


0.6 
0.4 
0.2 


0.6 
0.4 
0.2 




P(Y,K + )A 
cos(e K ) = o.05 



+ 



+ 



+ 



+ 




| l l l | l l l | l l l | 
COS(9 K ) = 0.45 



COS(9 K ) = 0.65 



I i i i I i i i I i i i I 
COS(9 K ) = 0.85 



p(y,K + )A 
cos(e K ) = o.i5 



+ 



+ 



+ 



+ 



COS(9 K ) = 0.35 



1 1 1 1 1 



+ 



1 1 1 1 1 



COS(9 K ) = 0.55 



COS(9 K ) = 0.75 



= '- t V'.i'Mt^f t ' M | , t ^ t t j t 



COS(e K ) = 0.95 



1.6 1.8 2 2.2 2.4 1.6 1.8 2 2.2 2.4 

W (GeV) W (GeV) 



0.4 
0.2 



0.4 
0.2 
^ 

.Q 

3 0.2 

a o 

~6 0.4 
"° 0.2 



0.4 
0.2 

1 



p(y,K + )^ 
- COS(6 K ) = o.os^^**^ 


■ ■ ■ I ■ ■ ■ I ■ ■ ■ I ■ ■ ■ I 

; p( y! K + )z°; 

- COS(9„) = 0.15 


_' ' ' M ' ' I ' ' ' M ' ' l_ 

COS(9 K ) = 0.25 


' ' I ' ' ' I ' ' ' I ' ' ' l_ 
" '^OS(9^)To*3T S * 4 * , * i:, ^ 


_i i i 1 ■ i i 1 i i i 1 ■ i i l_ 
" r COS(9 K ) = 0.45 " 


_i i i I i i i I i i i I i i i l_ 

- *^*^is(9^)Ta55* r< * r ""* - 


1 ■■■■■■■ 1 J 
COS(9 K ) = 0.65 

' ■ ■ ■ I ■ ■ ■ I ■ ■ ■ I ■ ■ . r 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 \'_ 

" COS(9 K ) = 0.75 

' i i i 1 i i i 1 i i i 1 i i i r 


COS(9 K ) = 0.85 


COS(9 K ) = 0.95 


... i ... i ... i ... i 



6 1.8 2 2.2 2.4 1.8 2 2.2 2.4 

W (GeV) W (GeV) 



FIG. 9: [color online] Comparison of the calculated and experimental differential cross sections for 
the p(7, K + )A and p(j, K + )TP reactions as a function of the 7p invariant mass (W) for various 
positive cosines of the K + cm. angle. The experimental data are taken from Ref. ij]. 



are well understood, this leads to the suggestion that a contribution to the background is 
still missing. It this context supplementing the high energy parts of the t— channel exchanges 



by Regge-trajectory exchange 64M68j may be an interesting option. Work is in progress to 
include the Regge-trajectory exchange in our model. 

For the backward angles (see Fig. 9), the data on the p^j, K + )A reaction show two peaks 
at W around 1.7 GeV and 1.9 GeV. Our calculations are able to reproduce the structure 
of both the peaks for all the angles except for the most backward one [COS(#^+) = 0.85] 
where the peak at W ~ 1.9 GeV is underestimated. We would like to stress that unlike 
Ref. {34)], we find no need of including an additional D Vi resonance (with mass « 1.9 GeV 
and width 0.316 GeV) to reproduce the data around 90°. The underestimation of the peak 
for W around 1.9 GeV for COS(6*ft-+) = 0.85 might be seen as an indication of the need to 
include such a resonance. However, with the inclusion of this resonance the cross section 
increases for other backward angles also around this value of W 



3s 



and the signs of the 



beam asymmetries come out to be opposite to what has been observed experimentally for 
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FIG. 10: [color online] Same as that shown in Fig. 9 but for various negative cosines of the K + 
cm. angle. The experimental data are taken from Ref. 3]. 



this W 

On the other hand, for the channel, the W dependence of the DCR (see Figs. 9 

and 10) is reproduced very well by our calculations for all the angles [corresponding to both 
negative and positive values of COS(9 K +)} with the exception of one very forward 9 K +. 

Polarization observables provide more sensitive tests of reaction models. The reason for 
this lies in the fact that these observables are generally very sensitive to the imaginary parts 
of the amplitudes which are governed by coupling to other channels via the optical theorem. 
Data on hyperon recoil polarization (Py, with Y being A or E°) have been reported by both 
SAPHIR and CLAS collaborations 2], yj- Py is related to interferences of the imaginary 
parts of the resonant amplitudes with real parts of other amplitudes including those of the 
back ground terms. In Fig. 11 we compare the results of our calculations Py with the 
corresponding data of the SAPHIR collaboration for ^(7, K + )A and p(7,i^ + )E° reactions. 
One notices that experimental Pa tend to be positive at backward angles, nearly zero at 
angles around zero and negative at forward angles. In contrast, F'yo data show nearly 
opposite trend. The CLAS Py data are used in the analyses of Refs. 34|, l38 |. 

Our calculations reproduce approximately the trends seen in the data. The opposite 
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FIG. 11: [color online] Comparison of our calculations with the experimental data (taken from 
Ref. [3]) for hyperon polarization for p(^,K + )h and p("y, K Jr )Y,° reactions as a function of the 
cosine of K + cm. angle for the photon energy bins (in GeV) as indicated in the figures. 

signs of the observed Py of two channels are nearly reproduced. The agreement with data 
is relatively better for E 1 below 1.6 GeV. The large positive back angle polarizations seen 
in the Pa data are reproduced for all the photon energy bins. Similarly the large negative 
experimental P^o at these angles are also nearly reproduced. At forward angles our model 
is relatively less successful in reproducing the Pa data at P 7 above 1.6 GeV. This again 
indicates perhaps some inadequacy of the background terms within our model. 

Beam asymmetry (S#) is the measure of the azimuthal anisotropy of a reaction yield 
relative to the linear polarization of the incoming photon. In Fig. 12 we compare the results 



of our calculations for this observable 
corresponding data taken from Ref. j|3 



or both p(7, K + )A and ^(7, K + )TP reactions with the 
6| which are available for nine energy bins for photon 
energies between 1.5 GeV to 2.3 GeV. Available experimental are positive for both the 
reactions. Our calculations reproduce the data with varying degrees of success. For the 
p(j, K + )A reaction the agreement with the data is relatively better for E 1 < 2.0 GeV. For 
larger E 1 the data are underestimated by our model. On the other hand, for the ^(7, K + )TP 
case the data are described better for E 1 > 1.9 GeV while at lower photon energies they are 



overestimate. In calculations reported in Ref. 



26] comparison with the data is shown for 
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FIG. 12: [color online] Beam asymmetry for the p(j, K + )A and p(7, K + )TP reactions as a function 
of the cosine of K + cm. angle for 15 photon energies. The experimental data are taken from 

n 

Ref. |5J. 

the K + A channel for one photon energy (E~ = 1.946 GeV) only. In the coupled-channels 



analysis of Ref. [34J, the agreement with the data for the p(j, K + )A reactions is of the 
same quality as that achieved by us - there too the data are better reproduced for E y < 
2.0 GeV while they are underpredicted at energies larger than this. However, we emphasize 
that such data for both K + A and K + TP channels have not been simultaneously analyzed 
with one parameter set in any other coupled-channels model. 



IV. SUMMARY AND CONCLUSIONS 

In this paper we investigated the photoproduction reactions p(7,_fT + )A and p( , y,K + )'E 
within a coupled-channels effective-Lagrangian approach which is based on the K-matrix 
method. Unitarity effects are correctly taken into account, since all important final channels 
(consisting of two-body systems irN, rjN, <pN, pN, 7V, KA, and KT,) are included in the K- 
matrix kernel. We build this kernel by using effective Lagrangians for the Born, w-channel, t- 
channel, and spin-^ and spin-| resonance contributions. Thus, the background contributions 
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are generated consistently and crossing symmetry is obeyed. The advantage of a full coupled- 
channel calculation is that it allows for the simultaneous calculation of observables for a large 
multitude of reactions with considerably fewer parameters than would be necessary if each 
reaction channel is fitted separately. More significantly, the implementation of unitarity 
ensures that the imaginary parts of the amplitudes are compatible with the cross sections 
for other channels. 

Our model provides reasonable description of the experimental data of the SAPHIR group 
on total and differential cross sections as well as on hyperon polarizations for photon energies 
ranging from threshold to up to ~ 3 GeV for both ^(7, K + )A and ^(7, K + )TP reactions. The 
beam asymmetry data of the SPring8 /LEPS group for the two reactions are also reasonably 
reproduced. 

An important point about our study is that the same parameter set was used in calculating 
observables for both the reactions. 

We made a detailed investigation of the contribution of various partial waves and showed 
that peaks seen in the total ^(7, A' + )A cross sections for photon energies around 1.5 GeV 
(invariant mass of 1.9 GeV) are largely due to coupled-channels effects rather than the 
contribution of a -D-wave resonance with mass around 1.9 GeV. A major part of cross 
sections of this reaction is generated via the background terms and resonances give prominent 
contributions only in the first peak region (for E 1 around 1.1 GeV). On the other hand, the 
total p(7, K + )TP cross sections are dominated by the P33 (1600) resonance in its peak region 
{E 1 around 1.5 GeV). 

The agreement between our calculations and the differential cross section data of the 
CLAS collaboration is not of the same quality as that seen in case of the SAPHIR data, 
particularly at extreme forward angles. At some backward angles the CLAS data show 
pronounced peaks for some photon energies but no peak is seen in the corresponding SAPHIR 
data or in our theoretical cross sections. Whether or not these peaks are an indication for 
an additional P or D wave resonance is a matter of debate and additional investigations. 
In this context, it is quite desirable to settle the issue of mutual inconsistency between the 
CLAS and SAPHIR data sets. 

Our work shows that it is indeed possible to fit meson photoproduction data of many 
channels simultaneously with a single parameter set within a coupled-channels model. Fur- 
ther improvement of our model is however, necessary for a better description of the data at 
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extreme forward angles in case of higher photon energies. 
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Appendix A: Effective Lagrangians 



We list here the effective Lagrangians for various vertices, p, k, pi and — q represent 
four momenta of the initial nucleon, final meson, final nucleon and photon, respectively. 
We assume that meson momenta are directed into the vertex, so that energy momentum 
conservation reads as p + k = pi — q. 

For the nucleon vertices the following couplings were used 



£>NNr) 
£NNcr 
C-NNp 

£nNuj 

£>NN<f> 



WNNn^N — : 75 * 



WNNri^N — j 75^7V 

X + 1 

-9NNa^N^fa^N 



N 



-gNNp^N [ipVp 1 + 
-gNN^N ( lp^ + 



c 



NN-y 



1 + T °-7^ 



C 



NN-yip 



gNN<t> 



2,171 j\r 

2,171 j\r 
K 



(Al) 



2m 



N 



2 2m^ 



'N 



The parameter x controls the admixture of pseudoscalar and pseudovector components in 
the corresponding Lagrangian. Its value is taken to be 0.5. This value was obtained in our 



previous study of photoproduction of associated strangeness [25j and has been held fixed 
in the study of all other reactions within our model. Nucleon spinors are depicted by \l/ 
and meson fields by <p. The magnetic moments are represented by k. Cnn^ generates the 



seagull or the contact term diagrams. We have followed the notations of Ref. 



69). 
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The Lagrangians for the meson vertices are 



Tfl-jr 



(A2) 



m p 

C PP1 = 2e(A»(d^jT <p p » - (d v A^ Pv r ^ 

+ (d»A^T Q <p Pw ) 
«-» 

<-> 

^tt,? = -ig fm r,(<Pr,d' x ip 7V )ip Pti 

777.^ 

The coupling constants entering into Eqs. (A.l) and (A. 2) together with baryon magnetic 
moments are listed in Table III. We have taken positive values for all primary coupling 
constants involving the nucleon. In particular, the sign of gNKA differs from its customary 



negative value 



70( However, we would like to stress that in a calculation like ours and also 



in many of those cited in Ref. [70| this sign is undetermined. Changing the sign of all the 
coupling constants involving a single A-field leaves the calculated observables invariant since 
it corresponds to a sign redefinition of this field. The magnitudes of the couplings are within 
the broad range specified in [tqI ] . 

For the Sn, Ssi,Pn and resonances the hadronic couplings are written as 

£ v nr 1/2 = -g<pNRVR\xir<p + (1 - X )j^^(d^)}^N + H.c, (A3) 
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TABLE III: Parameters summary table 





13.47 


gNNr) 


0.85 


gNNa 


10.0 


SNNp 


4.2 


SNNw 


3.0 


&NN<f> 


-0.0 


gEcrp 


33.0 


gSAp 


-27.0 
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0.0 




-0.12 


p- i 
opTv 'y 


-0.10 


OH'I ; 


-0.21 


gcj7^ 7 


-0.12 


guj7T 7 


0.32 


°r u / 


12.0 


gTVAK 


10.0 




14.5 


oiv iir\ 


-3.3 


Oiv Zj/V 


0.0 


g<^>KK 


-4.5 


gpKK 


-3.0 


g-n-KK* 


-3.26 


gjjKK* 


-3.2 


Sk*ko 7 


0.177 


gK*K±7 


-0.177 


Kp 


1.79 


K n 


-1.91 


K A 


-0.613 




0.79 


K S+ 


1.45 




-0.16 




-1.61 







where M = (rrif> ± m^), with upper sign for even parity and lower sign for odd parity 
resonance. The operator T is 75 and unity for even and odd parity resonances, respectively. 
For isovector mesons, ip in Eq. (A. 3) needs to be replaced by r ■ ip for isospin-^ resonances 
and by T • cp otherwise. 

The corresponding electromagnetic couplings are 

r 



C 



■yNR 1/2 



Am 



a^ N F^ + H.c, 



N 



(A4) 



where is the resonance spinor and F^ v = d^A v — d u A^. The operator T is 1 for the 
positive parity resonance and —275 for the negative parity one. 

For spin- 1 re sonances, we have used the gauge-invariant effective Lagrangians as discussed 



in Refs. 
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57 



5s 



71 



72] . We write here the vertex functions used by us in computation 
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involving these vertices. The resonance- nucleon-pion vertex function (e.g.) is given by 

~l a {q-p)-j><f\ [(1-X) + Xfi/M p ], 



■pet 



9i 



and the corresponding electromagnetic vertices are 



(A5) 



L R 3/2 ^N-y 



(92 + 2 9l 



q V - • q 



+ 



.91 



+ 



(A6) 



Here p is the four-momentum of the resonance and q is that of the meson. Index a belongs 
to the spin- 1 spinor and /x to photon. Interesting property of these vertices is that the 
product, p-Y — 0, where V defines the vertices on the left hand side of Eqs. (A. 5) and (A. 6). 
As a consequence, the spin-^ part of the corresponding propagator becomes redundant as 
its every term is proportional to either p^ or p v . Thus only spin-| part of this propagator 
gives rise to non-vanishing matrix elements. 
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